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in Hydrogen and Deuterium from 1.1 to 5.5 GeV 
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The differential cross sections for the 771 — -> %~p and the 7p — -> ir + n processes were measured at 
Jefferson Lab. The photon energies ranged from 1.1 to 5.5 GeV, corresponding to center-of-mass 
energies from 1.7 to 3.4 GeV. The pion center-of-mass angles varied from 50° to 110°. The tt~ and 
7r + photoproduction data both exhibit a global scaling behavior at high energies and high transverse 
momenta, consistent with the constituent counting rule prediction and the existing ir + data. The 
data suggest possible substructure of the scaling behavior, which might be oscillations around the 
scaling value. The data show an enhancement in the scaled cross section at center-of-mass energy 

near 2.2 GeV. The differential cross section ratios at high energies and high transverse 

momenta can be described by calculations based on one-hard-gluon-exchange diagrams. 

PACS numbers: 13.60.Le, 24.85. +p, 25.10.+S, 25.20.-x 



I. INTRODUCTION 

Quantum ChromoDynamics (QCD), a fundamental 
theory for describing the strong interaction, is not 
amenable to analytical solutions in the nonperturbative 
region. Some dynamical models must be developed. 
Meson-exchange models in terms of the nucleon-meson 
degrees of freedom describe nuclear physics data well at 
low energy, and perturbative QCD (pQCD) in terms of 
the quark-gluon degrees of freedom succeeds in explain- 
ing many measurements at high energy. But little is 
known about the transition between these two regions. 
Testing the constituent counting rule for the exclusive 
reactions is one way to study the transition of the de- 
grees of freedom. 

The constituent counting rule establishes a direct con- 



nection between the quark-gluon degrees of freedom and 
the energy dependence of the differential cross section for 
exclusive processes at fixed center-of-mass angles. This 
rule was first derived from simple dimensional count- 
ing y|,y,|2| and was later confirmed in a short-distance 
pQCD approach H. It is consistent with many exclusive 
measurements [ETH IS U El ■ However, there are still 
many puzzles. This rule begins to agree with experimen- 
tal data at photon energies as low as 1 GeV 0, whereas 
pQCD is not expected to be valid at such low energies. 
The hadron helicity conservation rule [TH ]. another out- 
come of the same short-distance pQCD framework, does 
not agree with data in the same energy and momen- 
tum transfer region 0, ^| . There are a few anomalies 
beyond the constituent counting rule in the extensively 
studied pp scattering process 0, 0, . Recently, the 
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parton orbital angular momentum has been found to play 
a non-negligible role in the exclusive reactions, which may 
explain hadron helicity nonconservation and other polar- 
ization measurements 0, . 

Single pion photoproduction, "fN — > ttN, is a rela- 
tively simple process for studying the strong interaction. 
It has larger cross sections at high energy than other ex- 
clusive channels due to its slower decrease with energy, 
i.e. da/dt ~ s~ 7 . One can also form the differential 

cross section ratio fffi''""" + P \ ■ The ratio is amenable 
to theoretical predictions since many factors may cancel 
out in leading order. 

This paper focuses on extracting the differential cross 
sections for the single charged pion photoproduction pro- 
cesses, 7p — ► TT + n and 771 — * 7r"p. This is one major 
goal of experiment E94-104 |l9j carried out at Jefferson 
Lab (JLab). The photon beam energy ranged from 1.1 
to 5.5 GeV. The pion center-of-mass angles varied from 
50° to 110°. The results at 90° have already been pub- 
lished [20j. which are also updated in this paper. The 
experiment E94-104 also measured ir~ photoproduction 
with a helium target, providing the first nuclear trans- 
parency data for this process |21j . 

This paper is organized as follows. Section II intro- 
duces the theoretical and experimental background for 
JLab experiment E94-104. Section III describes this ex- 
periment at JLab Hall A. Section IV presents the data 
analysis procedure to extract the differential cross sec- 
tions. Section V discusses the results. Section VI and 
VII are outlook and acknowledgments. 

II. THEORETICAL AND EXPERIMENTAL 
BACKGROUND 

The constituent counting rule is also called the dimen- 
sional scaling rule. It states that 

(da/dt) AB ^ C D ~ s- (n - 2) f(9 c . m .) (1) 

for an exclusive two-body reaction AB — > CD when 
s — > 00. Here s and t are the Mandelstam variables and n 
is the total number of elementary fields (quarks, leptons 
or photons) which carry finite fractions of particle mo- 
mentum. It states that at fixed center-of-mass angle B c . m . 
and large s, {da / dt) pp ^ P p ~ s~ 10 , (da / dt)^^^ ~ s~ 8 , 
(da/dt) yd-^pn ~ s 11 and (da/dt)^^^^^ ~ s ~ 7 . The 
constituent counting rule implies something of fundamen- 
tal importance: the quark has not only a mathematical 
existence, giving current algebra, Bjorken scaling and the 
hadron spectrum, but a dynamical existence as well |2|. 

The constituent counting rule was originally derived 
from simple dimensional counting by Brodsky and Far- 
rar p], and simultaneously by Matveev el al. Q in 1973. 
Brodsky and Farrar also examined the required condi- 
tions for the simple dimensional derivation: 
(a) the effective replacement of the composite hadron by 
constituents carrying finite fractions of the hadron mo- 
mentum, 



(b) the absence of any mass scale in the amplitude or 
binding corrections. 

They showed that both condition (a) and (b) are natural 
features of renormalizable field theories, with certain dy- 
namical assumptions concerning the nature of the wave 
function, the absence of infrared effects, and the accumu- 
lation of logarithms Q • 

Later in 1980, Lepage and Brodsky showed that the 
constituent counting rule can be reproduced within a 
short-distance pQCD approach , up to calculable pow- 
ers of the strong coupling constant. The energy depen- 
dence of the strong coupling constant is small at high 
energy, and some recent r decay data also su gges t the 
freezing of the coupling constant at low energy [22| . 

Another outcome of this approach is hadron helicity 
conservation fll| : 

h A + h B = h c + ho, (2) 

which leads to strong correlations between the final state 
helicities. The above results came from the calculation 
of an enormous number of connected tree diagrams for 
hard subprocesses without considering the parton orbital 
angular momentum, while the soft subprocesses, such as 
Landshoff diagrams [2j| were suppressed in leading order 
for example due to gluon radiation. 

The scaling behavior of the differential cross section 
predicted by the constituent counting rule can also be de- 
scribed by string theory [24) and other phenomenological 
models. For the deuteron photodisintegration process as 
an example, the Quark-Gluon String (QGS) model [2{| 
and the Hard Rescattering Mechanism (HRM) [2(j de- 
scribe fairly well both the energy dependence and the 
asymmetric angular distribution of the data 
The Reduced Nuclear Amplitudes (RNA) |29j and the 
Asymptotic Meson Exchange Calculation (AMEG) l3fJl 
can describe the scaling behavior at 9 c . m . = 90° [ll&IIJ 

On the experimental side, the constituent counting rule 
is consistent with data for many exclusive processes, such 
as ppelastic scattering [j| , hadron-hadron elastic scatter- 
ing |(J and deuteron photodisintegration 0, H, ■ The 
fitted power of - from the pp elastic scattering data with 
s > 15 GeV 2 and \t\ > 2.5 GeV 2 is equal to 9.7 ± 0.5 0, 
consistent with 10 as predicted by the constituent count- 
ing rule. Eight meson-baryon and two baryon-baryon 
exclusive reactions at 9 c . m . = 90° were measured at the 
AGS (the Alternate Gradient Synchrotron at BNL) with 
beam momenta of 5.9 GeV/c and 9.9 GeV/c. The fit- 
ted powers of - are also consistent with the constituent 
counting rule predictions, i.e. 8 for the meson-baryon re- 
actions and 10 for the baryon-baryon reactions, except for 
one reaction: n~p — ► 7r + A~ |(|. Deuteron photodisinte- 
gration is another process that exhibits scaling behavior 
of the differential cross sections 0, H, 13 ■ The onset of 
scaling in photon energy for deuteron photodisintegra- 
tion depends greatly on the center-of-mass angles. The 
corresponding threshold of the proton transverse momen- 
tum, Pt = J ^MdEySir^Oc.m.i is above about 1.1 GeV/c 
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for the proton angle between 30° and 150° 0. 

Despite the theoretical and experimental support for 
the constituent counting rule, there remain some puz- 
zles and anomalies. First of all, it is surprising to see 
the onset of scaling at transverse momentum as low 
as about 1.1 GeV/c 0, 0, such as in the photodisin- 
tegration data. The applicability of pQCD to exclu- 
sive processes remains controversial in the GeV region. 
The pQCD calculation fails to predict the magnitude of 
some fundamental quantities, such as the proton mag- 
netic form factor G P M with Q 2 up to 30 (GeV/c) 2 [3l[. 
Hadron helicity conservation, another consequence of 
pQCD (this statement is currently under debate |3^1. 
tends not to agree with polarization measurements, such 
as those from JLab for the photodisintegration process 
d(j,p)n up to 2.4 GeV [13] and neutral pion photo- 
production p("f,p)iT up to 4.1 GeV 0. Although 
contributions from nonzero parton orbital angular mo- 
mentum are power suppressed as shown by Lepage and 
Brodsky £| , they could break the hadron helicity con- 
servation rule |17|. Orbital angular momentum could 
also lead to asymptotic scaling of the proton form fac- 
tor ratio: F 2 (Q 2 )/F 1 (Q 2 ) ~ (log 2 Q 2 / A 2 ) / Q 2 with 0.2 
GeV< A <0.4 GeV based on an explicit pQCD calcula- 
tion [3 or F 2 (Q 2 )/F 1 (Q 2 ) ~ 1/^/Q 2 pl^l that agrees 
with the JLab proton form factor data 34] . A recent non- 
perturbative analysis |35| of the hadronic form factors 
based on light-front wave functions, and a model with 
an intrinsic (quark-like) structure and a meson cloud |3fil ] 
also describes the JLab proton form factor data [34[ well. 

Furthermore, several striking anomalies have been ob- 
served in pp scattering. One is the very large spin-spin 
correlation. The ratio of (da/dt)^ /(da/dt)n with spin 
normal to the scattering plane can reach 4 in pp elas- 
tic scattering at c . m . = 90° ^4|. Next is the oscillation 
of the differential cross section da/dt around the scal- 
ing value in pp elastic scattering [l5(- The third is the 
anomalous energy dependence of nuclear transparency of 
the A(p,2p) process fig . 

There exist different theoretical attempts to describe 
the anomalies. One example is the interference between 
two types of subprocesses, the short-distance hard sub- 
processes and the long-distance soft (Landshoff) subpro- 
cesses [33|. Another example is the interference between 
the pQCD background and two J = L = S=1,B = 2 
resonance structures associated with the strangeness and 
charm production thresholds [38j . 

Recently, more mechanisms were under discussion that 
can lead to the deviation from the scaling. For example, 
the generalized constituent counting rule [3JJ, I40l l4l| in- 
cluding the parton orbital angular momentum and the 
restricted locality of quark-hadron duality ^^] . 

Including the nonzero parton orbital angular momen- 
tum may change the expression of the constituent count- 
ing rule. Based on a hadronic light-cone wave function 
involving parton orbital angular momentum, which has 
been used to describe the JLab proton form factor data, 
a generalized constituent counting rule [39j can be de- 



rived for hard exclusive processes by counting the soft 
mass dimensions of scattering amplitudes. This general- 
ized constituent counting rule with parton orbital angu- 
lar momentum dependence can be also derived |40j in a 
non-perturbative method, which does not rely on pQCD. 
According to the generalized constituent counting rule, 
the fixed-angle scattering cross section behaves like 
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for the exclusive reaction A + B — ► C + D + where nn 
is the number of elementary fields in involved hadron 
H and Aa contains only angular variables. For par- 
ton orbital angular momentum I z h = 0, this is just the 
traditional constituent counting rule. As a result, the 
helicity-flip amplitudes for the pp — > pp process were pre- 
dicted to scale as s -9 / 2 with ^2 H \I z h\ — 1 or s~ 5 with 
y]fj \I z h\ = 2, while the helicity conserving amplitudes 
were known to scale as s~ 4 . The interference between 
amplitudes with different helicity changes offers a new 
mechanism to explain the spin-spin correlation and oscil- 
lation around the scaling value in pp scattering |4lj . This 
can also lead to the deviation from the traditional con- 
stituent counting rule for other exclusive processes, such 
as photoproduction of charged pions that is discussed in 
this paper. 

Therefore, a detailed investigation of the scaling be- 
havior may enable a test of the generalized counting rule, 
though a more rigorous test should come from the polar- 
ization measurements, which would allow the separation 
of amplitudes with different helicity changes. 

The deviation from the constituent counting rule for 
exclusive processes may also be due to the breakdown of 
the locality of quark-hadron duality [42^. Quark-hadron 
duality is an empirical property of the data discovered 
by Bloom and Gilman before the advent of QCD [43| . 
The production of resonances at lower energies and mo- 
mentum transfers averages smoothly around the scaling 
curve measured at large momentum transfers. The sum 
over resonances can be related to the scaling behavior as 
a result of destructive interference. This is rather local 
at high energy due to the high density of the overlapping 
resonances, which is called the locality of quark-hadron 
duality. But the local degeneracy may not be reached 
at energies of a few GeV, which leads to the restricted 
locality of quark-hadron duality. The restricted local- 
ity may cause oscillations around the scaling value above 
the resonance region when different partial waves are not 
canceled locally. 

The exclusive charged pion ratio of ^/jf[^p~^+n) 
can be estimated based on the one-hard-gluon-exchange 
Feynman diagrams |44l EE) . In Huang et a/.'s ap- 
proach 44], the helicity amplitude for the photoproduc- 
tion of a meson was assumed to factorize into the parton- 
level subprocess amplitude and the nucleon form factors. 
Evaluating the four Feynman diagrams in Figure Ogives 
the parton-level subprocess amplitude for pseudoscalar 
meson photoproduction. Due to isospin invariance, the 
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FIG. 1: One-hard-gluon-exchange Feynman diagrams for the 
parton-level subprocess jq — > Mq in the single meson photo- 
production *yN — > MN . 

form factors are divided out and the exclusive charged 
pion ratio takes on a simple form, 

da I dt(-yn — > 7r~ p) fued + se u \ 2 
da I 'dt(-yp — > 7T+ n) \ue u + sed ) 

where u and s are the Mandelstam variables and e q de- 
notes the charge of quark q. 

III. EXPERIMENT 

To study the transition from nucleon-meson degrees of 
freedom to quark-gluon degrees of freedom, it is essential 
to investigate the GeV region where the transition ap- 
pears to occur. While there were some measurements at 
SLAC for the TP — ► ir + n process at photon energies of 4, 5 
and 7.5 GeV [l(|, which exhibit a global scaling behavior 
expected by the constituent counting rule, there are no 
data beyond 2 GeV for 771 -> ir'p lM IItI |48 | . This ex- 
periment, JLab experiment E94-104 [ljj was proposed to 
measure the cross section for charged pion photoproduc- 
tion da / 'dt{^n — > 7r~ p) and da / dt(jp — > n + n) from 1.1 
to 5.5 GeV. In addition, the differential cross section ra- 
tio for charged pion photoproduction can be formed and 
compared to theoretical predictions. 

The experiment E94-104 was carried out in Jlab 
Hall A |5fJ, using the continuous electron beam at cur- 
rents around 30 /iA. The core of the Hall A equip- 
ment is a pair of nearly identical 4 GeV/c spectrome- 
ters capable of determining the momentum and angles 
of charged particles with high resolution. The schematic 
view of the setup for this experiment is shown in Fig- 
ure |21 The real bremsstrahlung photons were gener- 
ated by the electrons impinging on a copper radiator, 
located 72.6 cm upstream from the target. The foil with 
a thickness of 6.12% radiation length was used for the 
production data of E94-104. A liquid hydrogen target 
(LH2) was used as the proton target, while a liquid deu- 
terium target (LD 2 ) was used as an effective neutron tar- 
get. The outgoing pions and protons were detected by 
the two high resolution spectrometers (HRS) in Hall A. 
The vertical drift chambers (VDCs) recorded its track 
and the scintillator planes (S1/S2) provided timing in- 



formation and generated triggers. Aerogel Cerenkov de- 
tectors (A1/A2/AM) provided particle identification for 
positively charged particles, mainly pions and protons. 
Gas Cerenkov, preshower/shower detector and pion re- 
jector were used to discriminate negatively charged par- 
ticles, mainly electrons and pions. The left spectrome- 
ter was optimized to detect positively charged particles 
while the right one was optimized to detect negatively 
charged particles. However, both spectrometers had to 
contain detectors to identify both negatively and posi- 
tively charged particles, since there were a few reversed 
polarity kinematics. Based on two-body kinematics, the 
incident photon energies were reconstructed from the fi- 
nal states, i.e. the momentum and angle of the tt + in 
the singles measurements for the 7P — > 7r + n process, or 
momenta and angles of the n~ and p in the coincidence 
measurements for the 777. — > n~p process. 

The coincidence kinematics for the jn — > ir~p process 
are listed in Table [I] Normally, the negatively charged 
pions are detected by the spectrometer to the right of 
the beam line (viewed along the beam direction), and 
the protons detected by the left spectrometer. But a few 
kinematics require reversing the polarities of the spec- 
trometers, because the maximum momentum of the right 
spectrometer is only 3.16 GeV/c, while that of the left 
spectrometer is 4 GeV/c. The singles kinematics for the 
7P — * ir + n process are listed in Table [H] The positively 
charged pions were detected by the left spectrometer with 
positive polarity. The beam energies tabulated in Tabled 
and Table|n]were the nominal values used to set the spec- 
trometers, which may deviate from the measured ones 
by several MeV. The spectrometer momentum and angle 
settings were calculated by using a photon energy close to 
the beam energy, i.e. E e — 75 (MeV), where the multiple 
pion production processes were suppressed. 



IV. DATA ANALYSIS 

A. Overview 

The raw data from the data acquisition (DAQ) sys- 
tem were replayed or decoded by an event processing 
program, ESPACE (Event Scanning Program for Hall A 
Collaboration Experiments) using CERNLIB packages. 
The outputs were histograms and ntuples of physical vari- 
ables in the HBOOK format. The yield from the data 
was obtained by applying cuts on certain variables in the 
ntuples, such as trigger type, particle type, spectrome- 
ter acceptance and reconstructed photon energy. Next, 
the yield was normalized by beam charge and computer 
deadtime. To extract the differential cross sections, sim- 
ulations were carried out by using the modified MCEEP 
(Monte Carlo for (e, djp)) program j^jj written for JLab 
Hall A. The raw differential cross section (%)data was 
extracted by comparing the background subtracted yield 
from the data (Idata) with the yield from the Monte Carlo 
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FIG. 2: Schematic view of the experimental setup for E94-104. 



TABLE I: Spectrometer settings for coincidence kinematics. E e is the electron beam energy, 8 c . m . the pion center-of-mass 
angle, Pl (Pr) the central momentum for left (right) spectrometer with the sign indicating its polarity, and 9l {Or) is the 
central scattering angle for left (right) spectrometer. The Mandelstam variables \fs and —t are in the last two columns. 



l(GeV) 



Pl 
(GeV/c) 



9l 

o 



Pr 
(GeV/c) 



Q (GcV) (GeV/c 



1.173 
1.173 
1.173 



50.0 
70.0 
90.0 



+0.521 
+0.727 
+0.923 



60.17 
49.72 
39.75 



-0.953 
-0.838 
-0.706 



28.33 
41.45 
56.66 



1.71 
1.71 
1.71 



0.253 
0.467 
0.709 



1.721 
1.721 
1.721 



50.0 
70.0 
90.0 



+0.697 
+0.989 
+1.277 



58.32 
47.39 
37.37 



-1.433 
-1.238 
-1.015 



24.46 
36.02 
49.73 



1.99 
1.99 
1.99 



0.433 
0.798 
1.212 



1.875 
1.875 



50.0 
90.0 



+0.742 
+1.370 



57.79 
36.75 



-1.566 
-1.099 



23.64 
48.21 



2.06 
2.06 



0.484 
1.355 



2.558 
2.558 
2.558 
2.558 
2.558 
2.558 



50.0 
70.0 
90.0 
50.0 
70.0 
90.0 



+0.913 
+1.322 
+1.740 
-2.108 
-1.794 
-1.438 



55.67 
44.37 
34.45 
20.96 
31.02 
43.18 



-2.108 
-1.794 
-1.438 
+0.913 
+1.322 
+1.740 



20.96 
31.02 
43.18 
55.67 
44.37 
34.45 



2.35 
2.35 
2.35 
2.35 
2.35 
2.35 



0.696 
1.282 
1.948 
0.696 
1.282 
1.948 



3.395 
3.395 
3.395 
3.395 
3.395 



50.0 
70.0 
90.0 
100.0 
110.0 



+1.113 
+1.642 
+2.195 
+2.466 
+2.725 



53.21 
41.74 
32.01 
27.69 
23.65 



-2.799 
-2.363 
-1.866 
-1.614 
-1.369 



18.57 
27.56 
38.57 
45.24 
53.01 



2.67 
2.67 
2.67 
2.67 
2.67 



0.971 
1.789 
2.718 
3.190 
3.648 



4.232 
4.232 
4.232 



50.0 
70.0 
90.0 



-3.489 
+1.949 
+2.638 



16.84 
39.51 
30.01 



+1.300 
-2.929 
-2.291 



51.04 
25.05 
35.18 



2.95 
2.95 
2.95 



1.248 
2.299 
3.494 



5.618 
5.618 



70.0 
90.0 



-3.863 
+3.359 



22.08 
27.38 



+2.442 
-2.990 



36.48 
31.11 



3.36 
3.36 



3.148 
4.785 
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TABLE II: Spectrometer settings for singles kinematics. E s is the electron beam energy, 9 c . m . the pion center-of-mass angle, 
Pl the central momentum for left spectrometer with the sign indicating its polarity, and 8l is the central scattering angle for 
left spectrometer. The Mandelstam variables y/s and —t are in the last two columns. 
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A 

"cm. 


r> 

rh 


f\ 

UL 


J* 


-t 


( t~y ~\ t\ 
(Gev) 
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(GeV/c) 
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(GeV) (GeV/c) 2 


1.173 


70.0 


+0.838 


41.45 


1.71 


0.467 


1.173 


90.0 


+0.706 


56.66 


1.71 


0.709 


1.721 


50.0 


+1.433 


24.46 


1.99 


0.433 


1.721 


70.0 


+1.238 


36.02 


1.99 


0.798 


1.721 


90.0 


+1.015 


49.73 


1.99 


1.212 


1.875 


50.0 


+1.566 


23.64 


2.06 


0.484 


1.875 


90.0 


+1.099 


48.21 


2.06 


1.355 


2.558 


50.0 


+2.108 


20.96 


2.35 


0.696 


2.558 


70.0 


+1.794 


31.02 


2.35 


1.282 


2.558 


90.0 


+1.438 


43.18 


2.35 


1.948 
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2.67 


0.971 


3.395 


70.0 


+2.363 


27.56 


2.67 


1.789 


3.395 


90.0 


+1.866 


38.57 


2.67 


2.718 


3.395 


100.0 


+1.614 


45.24 


2.67 


3.190 


3.395 


110.0 


+1.369 


53.01 


2.67 


3.648 


4.232 


70.0 


+2.929 


25.05 


2.95 


2.299 


4.232 


90.0 


+2.291 


35.18 


2.95 


3.494 


4.232 


100.0 


+1.967 


41.36 


2.95 


4.101 


5.618 


90.0 


+2.990 


31.11 


3.36 


4.785 


5.618 


100.0 


+2.547 


36.69 


3.36 


5.615 



target along the direction perpendicular to the spectrom- 
eters. The 8 is the relative deviation from the central 
momentum setting. 

Six two-dimensional boundaries are defined for each 
spectrometer, out of any two combinations of the four 
acceptance variables, 9 tg , 4>t g , ytg and S. Each bound- 
ary is a polygon defined in a two-dimension plot of the 
data. For each event, the magnitude of the distance to 
the boundary is normalized based on the maximal length. 
The sign of the distance is positive for the events inside 
the polygon. The R-function for a single spectrometer, 
useful for singles jp — ► tt + u data, is defined to be the 
minimal distance to the six boundaries, while that for two 
spectrometers, useful for coincidence 771 — ■> n~p data, is 
defined by twelve two-dimension boundaries. 



C. Particle Identification Analysis 

For the detection of positively charged particles, the 
protons need to be selected for the coincidence "fn — > ir~p 
process and the pions need to be selected for the sin- 
gles 7/) — > ir + n process. The aerogel Cerenkov detectors 
A1/A2/AM (see Figure |2J) were utilized to identify pro- 
tons and pions. The Al and A2 detectors were used for 
normal polarity data, while AM was used for reversed 
polarity data. The particle identification with aerogel 
detectors was consistent with other methods available 
at low momentum, for example measuring the Time-Of- 
Flight of the particle or the energy deposited in the scin- 



simulation (Y mc ): 

(dcr _ ,dcr. Vdata 

l di )data - ( dP mc *1^7- (5j 

Finally to extract the physical differential cross section, 
corrections such as the nuclear transparency of deuterium 
for the final state interaction, the detection efficiency and 
nuclear absorption in the detection materials were ap- 
plied to the raw differential cross section. 



B. Acceptance Analysis 

The R-function is defined to be the minimal distance 
to the acceptance boundary in terms of several two- 
dimensional polygons. It helps to select events in the 
central region of the spectrometer acceptance in a sys- 
tematic and efficient way, where the optics matrix ele- 
ments are well tuned. This method was originally devel- 
oped in the E89-044 data analysis |5^[. The version |53| 
refined in the E91-011 data analysis was used to analyze 
the E94-104 data and will be discussed below. 

The R-function is generated for each event in both data 
analysis and simulation to optimize the cuts on different 
acceptance variables, i.e. 8 tg , 4>tg, ytg and 5. The 9 tg 
is the deviation of the out-of-plane angle from the spec- 
trometer central setting of zero. The <f>t g is the deviation 
of the in-plane angle from the central setting of scatter- 
ing angle. The ytg is the reaction vertex position in the 
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tillators. 

For the detection of negatively charged particles, the 
pions need to be identified from the electron background 
for the coincidence "fn — > n~p process. It was realized by 
using the combination of the gas Cerenkov and shower- 
type detector, i.e. the gas Cerenkov and total shower 
(preshower/shower) detector for normal polarity kine- 
matics, the gas Cerenkov and pion rejector detector for 
reversed polarity kinematics. A one-dimensional cut was 
used for gas Cerenkov detector, while a two-dimensional 
graphic cut was defined for the shower-type detector. 

D. Background Subtraction 

The production data contained various kinds of back- 
ground, such as those from the electroproduction process 
and those from the end caps of the target. Therefore each 
complete kinematics consisted of four different configu- 
rations for data taking, 

(1) Radiator in, production target, 

(2) Radiator in, background target, 

(3) Radiator out, production target, 

(4) Radiator out, background target. 

For the coincidence 771 — * n~p process, the LD2 target 
was used as the production target and the LH 2 target for 
the background subtraction. For the singles jp — > Tr + n 
process, the LH 2 target was used as the production target 
and the dummy target for the background subtraction. 

The backgrounds were subtracted from the coincidence 
7r~ or singles ir + production yield according to 

71" : (51n,LD 2 — ^in,LH 2 ) ~ f (E 1 )(Y out: LD 2 ~ ^out,LH 2 ) 

7T : (^in,LH 2 ^in,Dum) /(-^7)(^out,LH 2 ^out,Dum )(6) 

where the in or out represents using or removing the ra- 
diator during the data taking. The factor f{E~ f ) is less 
than unity due to the interaction between the radiator 
and the electron beam. For most kinematics, the yield 
without the radiator was about one third of that with 
the radiator. The yield without the production target, 
especially for coincidence cases, was much smaller than 
the production yield. 

E. Monte Carlo Simulation 

The Monte Carlo simulation was performed by using 
MCEEP |51|. a computer program designed for coinci- 
dence (e, e'X) experiments in Hall A. The MCEEP pro- 
gram employs a uniform random sampling method to 
populate the experimental acceptance. An event is de- 
fined as one combination of variables that completely 
specifies the reaction in the laboratory. The cross section 
is considered as the weight of the event. It was modified 
for the coincidence 71(7, 7r~p) and singles p(j, ir + )n pro- 
cesses. The program was also modified to generate the 
R-function to define the acceptance cuts. 



The momentum distribution of the neutron inside the 
deuterium target was considered in calculating the kine- 
matics and cross section. The bremsstrahlung photon 
energy was randomly generated and the bremsstrahlung 
photon yield spectrum was calculated using thin-radiator 
calculation based on reference [54[, which was later cor- 
rected by the thick-radiator calculation [54l I5BI |5(j . The 
pion survival factor was also included in the cross section 
calculation based on the average flight length and was 
later corrected by muon contamination. 

The differential cross section % for n + photoproduc- 
tion at fixed center-of-mass angle was assumed to be 

dt = (I + cos0 c . m .) 4 (l - cos0 c . m .) 5 ' [ ~> ' U 

where the angular distribution was fitted to SLAC data in 
the several GeV region ^(J- The factor of 0.69 and so = 
10.263 came from these SLAC data at 5 GeV. The s~ 7 
dependence was consistent with the constituent counting 
rule. For 7r~ photoproduction, there was another con- 
stant to account for the deviation from tt + photoproduc- 
tion, but the overall normalization factor did not affect 
the final cross section extraction. The pion survival fac- 
tor Survive was calculated by F survivc = exp(- £^ c ) , 
with the average flight length L of 25 m and pion mean 
lifetime of r w = 2.60 x 10~ 8 s. 

The distributions of different variables from the data 
were compared with those from the Monte Carlo sim- 
ulation, including the acceptance variables, the recon- 
structed photon energy, the reconstructed initial momen- 
tum distribution of neutron inside the deuterium target, 
and the reconstructed center-of-mass angle, as shown in 
Figure |3 to Figure The results from data are symbol- 
ized by solid circles for coincidence measurements or by 
solid squares for singles measurements, while those from 
simulation are plotted as lines. Each simulated spectrum 
was multiplied by an arbitrary normalization factor for 
the comparison. This normalization factor is directly re- 
lated to the cross section extraction, as shown in Equa- 
tion [SJ Except for the photon energy comparisons, only 
one typical kinematic setting is shown for other cases, 
i.e. the one at E e = 4236 MeV and 9 c . m . = 90° for the 
coincidence measurements, or at E e — 4236 MeV and 
6c.m. = 90° for the singles measurements. 

The overall agreement between data and simulation is 
good, but there is some discrepancy for some cases. Since 
only reaction of interest was considered in the simulation, 
some discrepancy may be due to other physical reactions 
such as the multi-pion production. But the contribution 
from the multi-pion production is negligible with the cut 
on photon energy, as shown by the shaded area in Fig- 
ure |5l and Figure The uncertainties associated with 
various measurements such as the acceptance will affect 
the agreement too, which was included into the system- 
atic uncertainties. 
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FIG. 3: Comparison of acceptance variables from left and 
right spectrometers between data and simulation for coinci- 
dence measurements. 



F. Corrections and Uncertainties 

Table imi lists the systmatic uncertainties for kinemat- 
ics at 9 c . m , — 90°. The kinematics at other angles have 
similar breakdowns. 



TABLE III: The systematic uncertainties for coincidence and 
singles measurements at c . m . = 90°. 



Kinematics 


Coincidence (%) 


Singles (%) 


Photon Yield 


3.0 


3.0 




1.0 


1.0 


Target Density 


1.0 


1.0 


Tracking 


3.0 


2.0 


Scintillator/Trigger 


1.3 


1.0 


PID 


2.2 


3.0 


Muon Contamination 


3.0 


3.0 


Transparency 


5.0 




Nuclear Absorption 


4.2 


3.0 


Trial Cross Section 


2.2 


2.2 


Deadtime 


2.0 


2.0 


Acceptance 


3.0 


3.0 


Energy Loss 


2.0 


2.0 


Random Coincidence 


1.0 




Momentum Distribution 


1.0 




Beam Charge 


1.0 


1.0 


Total 


10.1 


8.1 



The thin-radiator calculation embedded in the simu- 
lation to calculate the photon flux was corrected for the 




-5 -2.5 2.5 5 

8 tg (%) 



-50 -25 25 50 
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FIG. 4: Comparison of acceptance variables from left spec- 
trometers between data and simulation for singles measure- 
ments. 



radiator thickness, by using a thick-radiator code |55j 
based on references 0; , which considers the energy 
loss of electrons in the radiator and is expected to be 
accurate to a level of 3% for the energies and radiator 
thickness used in this experiment. For the 6.12% copper 
radiator, the thin-radiator calculation overestimates the 
yield by 11% to 20%. By using this code 55], the f{E 1 ) 
can be also calculated, which will affect the subtraction 
of electron induced background. 

The decrease of target density due to beam induced 
local boiling was corrected. According to a comparison of 
the normalized yields with different beam currents, this 
correction was proportional to the average beam current 
and can be parameterized as 0.072* J (%) for LD2 target 
and 0.048 * I (%) for the LH 2 target with I the beam 
current in the unit of fiA. The uncertainty was on the 
level of 1%. 

The single wire efficiency of the VDCs, with samples 
defined by two neighboring wires, was very close to 100%. 
But multi-track events, dominated by two-track events, 
may still cause inefficiency in tracking. The corrections 
were applied based on the ratios of multi-track to single- 
track events, which were less than 2% for most kinemat- 
ics. The additional inefficiency associated with the track 
reconstruction algorithm was included in the uncertainty 
of 3%. 

The scintillator/trigger efficiency was obtained by 
checking the trigger for those selected events with good 
signals in VDCs and PID detectors. Special data were 
take to measure the efficiency, which was averaged to be 
98.8±1.1% for the right spectrometer and 98.8±0.7% for 
the left spectrometer. 
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FIG. 5: Comparison of reconstructed photon energy be- 
tween data and simulation for coincidence measurements at 
Oc.m, = 90°. The results from data are plotted as symbols, 
while those from simulation are plotted as lines. The electron 
beam energies are 1173.3, 1723.4, 2561.5, 3400.0, 4236.4 and 
5614.4 MeV. The comparison at beam energy 1876.9 MeV 
(not shown here) is very similar to that at 1723.4 MeV. The 
shaded events were chosen to extract the differential cross 
section. 
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FIG. 6: Comparison of reconstructed photon energy between 
data and simulation for singles measurements at 9 c . m . = 90°. 
The results from data are plotted as symbols, while those from 
simulation are plotted as lines. The electron beam energies 
are 1173.3, 1723.4, 2561.5, 3400.0, 4236.4 and 5614.4 MeV. 
The comparison at beam energy 1876.9 MeV (not shown here) 
is very similar to that at 1723.4 MeV. The shaded events were 
chosen to extract the differential cross section. 



One need to identify protons and negative pions for 
coincidence ir~ photoproduction, and identify only pos- 
itively charged pions for singles 7r + photoproduction. 
There is no correction on proton identification for coin- 
cidence ir~ photoproduction due to the very high p/tt + 
ratio ( > 100 ). The uncertainty due to 5-electrons (or 
knock-on electrons) was estimated to be 1%. For pion 
selection from coincidence tt~ photoproduction and sin- 
gles 7r + photoproduction, two PID detectors were used. 
The correction and uncertainty were estimated based on 
the performance study of each detector by defining good 
event samples based on the other one. The corrections 



depended on the particle momentum and the signal to 
noise ratio, and therefore varied by kinematics. 

The pion decay loss was considered as the survival fac- 
tor in the simulation by using the pion flight length and 
pion lifetime. The survival factor ranges from 53% to 
89% for different pion momentum. However, since some 
of the muons from pion decay may still fall into the accep- 
tance and be misidentified as pions, the calculation above 
may underestimate the effective pion survival factor and 
should be corrected. Based on the estimation with the 
modified simulation program SIMC ]^% . the correction 
depended on the particle momentum and ranged from 
-7% to -4% approximately [5^. The uncertainty was 
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FIG. 7: Comparison of reconstructed momentum distribution 
of the neutron in the deuterium target between data and sim- 
ulation for coincidence measurements. 
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FIG. 8: Comparison of reconstructed center-of-mass angle 
between data and simulation for coincidence measurements. 
The nominal center-of-mass angle is 90 degrees. 
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FIG. 9: Comparison of reconstructed center-of-mass angle be- 
tween data and simulation for singles measurements. The 
nominal center-of-mass angle is 90 degrees. 



estimated to be 3% by checking the dependence of the 
correction on the acceptance cuts. 

Nuclear effects must be considered to obtain the 
cross section for ^/n — > n~p from the measurement of 
d( / y,ir~p)p. The cross section may be reduced due to 
the final state interactions with the spectator proton in- 
side the deuteron. The nuclear transparency was defined 
to be the ratio of the reduced cross section to the raw 
cross section without any final state interactions. The 
measured nuclear transparencies for c?(e, e'p) quasi-elastic 
scattering 59], show little Q 2 dependence above Q 2 ~ 2 
(GeV/c) 2 ) and agree well with a Glauber calculation [r30| . 
The fitted value of 0.904±0.013 was used to deduce the 
nuclear transparency for the experiment E94-104, based 
on the Glauber formulation. The nuclear transparency 
for d('y,ir~p)p was scaled from that of d(e,e'p)n by re- 
placing the total pn scattering cross section with the to- 
tal pp and tt~p scattering cross sections. The systematic 
uncertainties in the nuclear transparency were estimated 
to be 5% to account for the uncertainties in the nuclear 
transparency measurement for c?(e, e'p)n and those in the 
effective pn, pp and n~p scattering cross sections. 

The produced particles, pions and protons in the coin- 
cidence measurements and pions in the singles measure- 
ments, had to go through various materials in the target 
and spectrometers before being detected. The event loss 
in the material is called nuclear absorption here. The ma- 
jor sources of nuclear absorption for high energy protons 
are listed in Table I1VI The absorption was calculated 
based on the thickness and effective absorption length of 
the material in the flight path of the produced particles. 
The effective absorption length A was estimated from the 
nuclear collision length At and nuclear interaction length 
A/ [61j as 2AtA//(At + A/) by assuming that half of the 
elastic and quasi-elastic scattering contribute to the ab- 
sorption. Later, the nuclear absorption was adjusted due 
to different flight lengths in the target, and different effec- 
tive absorption lengths for pions and protons at various 
momenta. The flight length can be calculated from the 
scattering angle and the geometry of the target. The en- 
ergy dependence of the effective absorption length was 
obtained from the cross section data in Reference |6lj . 
The uncertainty for each produced hadron was estimated 
to be 3%. 

The final cross section, extracted by comparing data 
and simulation, depended on the angular distribution and 
energy dependence of the trial cross section as the input 
of the simulation. Instead of searching for the exact form 
of the angular distribution and energy dependence of the 
actual cross section, the cross section fitted to SLAC data 
at high energy [lfj ] was used and its deviation from the 
actual cross section was considered as the systematic un- 
certainty. The trial angular distribution in the simula- 
tion may not be the same as the real case. The result- 
ing systematic uncertainty was estimated by checking the 
change of the final cross section from the flat angular 
distribution. The changes of the final results were very 
small due to the small acceptance of the Hall A spec- 
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TABLE IV: Major nuclear absorption in the target and spec- 
trometer for high energy protons. 



Material 


Thickness Density 


A 


Absorp. 




(cm) 


(g/cm 3 ) (g/cm 2 ) 


(%) 


15cm TjDo ci9Ki 


5.5 


0.1670 


49.8 


1.84 


15cm LH 2 (22K) 


5.5 


0.0723 


46.8 


0.85 


air 


300. 


1.21e-3 


73.4 


0.49 


SI (Polystyrene) 


0.5 


1.032 


68.3 


0.76 


S2 (Polystyrene) 


0.5 


1.032 


68.3 


0.76 


Al (Aerogel) 


9.0 


0.060 


75.5 


0.72 


A2 (Aerogel) 


5.0 


0.220 


75.5 


1.46 


AM (Aerogel) 


9.0 


0.100 


75.5 


1.19 


Gas Cerenkov (CO2) 


150. 


1.98e-3 


73.6 


0.40 



trometers. In the comparison, the mean center-of-mass 
angles were determined from the data, which deviated 
from those determined from the simulation by 0.2° at 
most. The systematic uncertainties due to the angular 
distribution were estimated to be 1%, 2% and 3% for 
kinematics at 9 c . m . = 90°, 100°, at 6 c . m . = 70°, 110°, and 
at c .m. = 50° respectively. The trial cross section used 
in the simulation had as -7 energy dependence, which 
was suggested by the SLAC data at high energy [l(| and 
was predicted by the constituent counting rule. The ac- 
tual energy dependences for both coincidence and singles 
kinematics are shown in Section The data, especially 
at low energy, do not have the s~ 7 energy dependence. 
The resulting systematic uncertainty was assigned to be 
2%, by checking the change of the final cross section from 
the flat energy dependence. 

The computer deadtimc was calculated by taking the 
ratio of missed triggers in the data stream from DAQ 
to the input triggers from scalers. It was less than 20% 
for nearly all the data and was corrected run-by-run (not 
listed here). The uncertainty was estimated to be around 
10% of the correction. The electronics deadtime was less 
than 0.5% for the majority of E94-104 data based on 
measurements using test pulses. There were also other 
systematics uncertainties for example due to acceptance, 
energy loss, random coincidence subtraction and beam 
charge. 



G. Reversed Polarity Data 

The polarities of the spectrometers were reversed for a 
few kinematics during the experiment. The particle iden- 
tification was optimized for data acquisition with nor- 
mal polarities. The aerogel detector (AM) in the right 
spectrometer was not as good as the combination of two 
aerogel detectors (Al and A2) in the left spectrometer 
in identifying protons. But since the proton signals are 
very clean with very low pion background, this hardly 
affected the results. The gas Cerenkov detector and pion 
rejector in the left spectrometer did not perform as well 



as the gas Cerenkov detector and preshower/shower de- 
tector in the right spectrometer in identifying pions. A 
tighter cut on momentum (—4% < S < 0%) was applied 
in the data analysis to avoid using the bad PMTs in the 
left gas Cerenkov detector. The pion rejector was only 
used to estimate the corrections. The reconstructed pho- 
ton energy spectrum from the data agreed with that from 
the simulation, as shown in Figure ITU1 For a consistency 
check, there were also some data recorded with both nor- 
mal polarities and reversed polarities. The differences in 
yields were within 5%, smaller than the systematic un- 
certainties (on the level of 10%). 
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FIG. 10: Comparison of reconstructed photon energy between 
data and simulation for coincidence 7r~ photoproduction with 
reversed polarities. The plot on the left represents the kine- 
matics at E e = 4236.4 MeV and 6 c . m . = 50°, while the plot 
on the right represents the kinematics at E e = 5614.4 MeV 
and 9cm. = 70°. 



V. RESULTS AND DISCUSSION 

The differential cross sections da/dfl and s 7 da/dt ex- 
tracted from JLab experiment E94-104 are shown in Ta- 
bleland Table EH for different beam energies and pion 
center-of-mass angles. The published results at 90° [2(j 
are also updated here. The angular distributions for all 
the energies are plotted in Figure ITT1 and Figure IT2l Also 
plotted are the SLAC tt+ data at E 1 = 4, 5, 7.5 GeV pH • 
Both the 7r + and 7r~ data at comparable energies, i.e. 
Ey = 4.2, 5.5 GeV, are consistent with the fit of the high 
energy SLAC data. 

It is worth mentioning that the pQCD calculations in- 
volving gluon self-coupling |62j cannot reproduce the an- 
gular distribution of the SLAC data, especially at the 
backward angles. This discrepancy may be due to the 
relatively low values of s, \t\ and u. The main contam- 
ination of their leading-twist predictions came from the 
i-channel meson resonances at forward angles, and from 
the u-channel baryon resonances at backward angles. 
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TABLE V: Differential cross sections for the — » n p process followed by the statistical and systematic uncertainties. 
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TABLE VI: Differential cross sections for the 7p — > ir + n process followed by the statistical and systematic uncertainties. 
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A. Comparison with the World Data 

As shown in Figure ^5] and Figure El the data from 
JLab experiment E94-104 (in solid circles) extended the 



singlepion photoproduction measurements at several 
GeV [lCj, Ha, H3, H3 > by spanning the resonance region 
and the scaling region. The differential cross sections of 
the jn — > Tr~p process with greater than 2.2 GeV 
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FIG. 11: Angular distributions from JLab E94-104 for the 
7?i — > ir~p process, as well as those from the SLAC data |To|| 
for the 7p — > -K + n process at photon energy of 4 GeV (open 
squares), 5 GeV (open circles), and 7.5 GeV (open triangles). 
The curve in each panel is the empirical fit of SLAC data: 
0.828e7- (1 - cos0 c .m.)" 5 • (1 + cos(9 c . m .)~ 4 . 
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were measured for the first time. The uncertainty in ^fs 
due to the 100 MeV photon energy window ranges from 
0.05 to 0.03 GeV as the beam energy increases 

The data agree within uncertainties with the world 
data in the overlapping energy region, except with the 
Besch et al. data |47j (open triangles). The Besch et al. 
data from Bonn suggest a very sharp peak in the scaled 
cross section for the jn — > ir~p process with t/s around 
2.0 GeV. Our data confirm the scaled cross section en- 
hancement around that region, but the peak is much less 
pronounced. We do not know the origin of this discrep- 
ancy exactly, though our momentum resolution (0.02%) 
is much better than that of Besch et al. (4%). The broad 
structure suggested by our data is seen in 7r° photopro- 
duction channel as well |46j | . Similar broad structure was 
also seen in the ir + and ir° channels from the preliminary 
JLab CLAS (CEBAF Large Acceptance Spectrometer) 
results l63T. 



y n — > ji p 
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FIG. 13: (Color online). Scaled differential cross section s 7 4f 
versus center-of-mass energy y^i for the yn — > tt~ p p rocess 
from JLab E94-104 and previous world data |4^ . l4?L l4^ | . The 
open triangles in the upper panel are averaged from the Besch 
et al. data [47| at 9 c . m . = 85° and 95°, while open triangles 
in the middle panel are averaged from those at 8 c . m . ~ 65° 
and 75°. The Fujii et al. and Schemer et al. data |48| in the 
middle panel were taken at 8 c . m . = 75°. 



FIG. 12: Angular distributions from JLab E94-104 for the 
7P — > 7r + n process, as well as those from the SLAC data [TH ] 
for the 7p — > n + n process at photon energy of 4 GeV (open 
squares), 5 GeV (open circles), and 7.5 GeV (open triangles). 
The curve in each panel is the empirical fit of SLAC data: 
0.828e7 ■ (1 - cos(9 c . m .)~ 5 • (1 + cos(9 c . m .)- 4 . 



B. Scaling at High Energy 

Based on the constituent counting rule, the differen- 
tial cross section da/dt at a fixed center-of-mass angle 
for the jn — > ir~p and 7p — > n + n processes is predicted 
to scale as s -7 . For both ir~ and tt + photoproduction 
processes, the data with y/s > 2.7 GeV at c . m . = 90° 
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FIG. 14: (Color online). Scaled differential cross section s 7 ^ 
versus center-of-mass energy yfs for the 7p — > 7r + n process 
from JLab E94-104 and previous world data [13, EES- 
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FIG. 15: (Color online). Scaled amplitude P^|M| versus 
center-of-mass energy y/s and transverse momentum Pt for 
the 7n — > n~p process. All the data points came from JLab 
E94-104. The arrows indicate the position below which the 
data do not scale for either 7r + or ir~ photoproduction. 



and y/s > 3.0 GeV at 8 c . m . — 70° indicate the scaling 
behavior of the differential cross section predicted by the 
constituent counting rule. The fitted power of - from 
the three data points in this region was 6.9±0.2 for the 
jn — * Tt~p process and 7.1±0.2 for the "fp — > ir + p pro- 
cess, consistent with the prediction of 7. This may have 
some theoretical implications, for example the validity of 
quark-gluon degrees of freedom and the freezing of the 
running strong coupling constant at several GeV. 

There is no sign of s~ 7 scaling for the data at 9 c . m . = 
50° up to center-of-mass energy of 3.0 GeV for the ir~ 
case and 3.9 GeV for the tt + case. This is not surpris- 
ing since the deuteron photodisintegration data 0, H 
at forward angles do not scale at as low energies as 
those at 90°. The deuteron photodisintegration data at 
Qc.m. — 53° seem to scale when the photon energy is 
greater than 3 GeV while the data at 9 c . m . = 90° scale 
when the photon energy is greater than 1 GeV. The cor- 
responding center-of-mass energies are 3.8 GeV and 2.7 
GeV respectively. 

The scaled invariant amplitude P^|M| is plotted in 
Figure ^] and Figure EH against center-of-mass en- 
ergy y/s and transverse momentum Pt, similar to what 
was performed for neutral pion photoproduction on the 
deuteron [64| • The invariant amplitude M was calculated 
from the differential cross section by using 



y p -> 7i n 



\M\ 41, ,„ : ; uJtt^^N^ttN) 



(8) 



and the transverse momentum was calculated by using 
Pt = | Pit I sm $7r from the pion momentum p^ and scat- 
tering angle 9^. The scaled amplitude P^|M| is plotted 
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FIG. 16: (Color online). Scaled amplitude P^|M| versus 
center-of-mass energy y/s and transverse momentum Pt for 
the 7p — > iv + n process. The data points in solid circles came 
from JLab E94-104. while those in open triangles are SLAC 
data [Tc|| . The arrows indicate the position below which the 
data do not scale for either 7r + or ir~ photoproduction. 



here for the jN — > ir^N process, while P^l-M] was plot- 
ted in reference [64| for the — ■> ir Q d process. The 
scaling power of Pt can be estimated by dimensional 
counting. Since the arrows in different colors (indicat- 
ing the possible onset of scaling for different angles) are 
closer in terms of transverse momentum than in terms 
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of the center-of-mass energy, the transverse momentum 
may be a better choice to describe the scaling onset than 
the center-of-mass energy, which was also stated previ- 
ously 0, 12^, [65|. The photoproduction data seem to 
reach the scaling region when the transverse momentum 
is around 1.2 GeV/c. As a comparison, the deuteron pho- 
todisintegration data start to exhibit scaling when proton 
transverse momentum ranges from 1.0 to 1.3 GeV/c at 
proton center-of-mass angle between 30° to 150°, except 
that the proton transverse momentum threshold ranges 
from 0.6 to 0.7 GeV/c at center-of-mass angle of 45°, 
135° and 145° [jj. Another interesting observation is 
that the scaled amplitude has much less dependence on 
center-of-mass angle than do the scaled cross sections. 
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FIG. 17: (Color online). Scaled differential cross section 
s? 1? versus center-of-mass energy y/s for the yn — > ir~p 
process (upper panel) and 7p — > ir + n process (lower panel) 
at 9 c .m. = 90° from JLab E94-104 and previous world 
data 0, ■ The error bars in the insets include 

only point-to-point uncertainties to highlight the possible sub- 
structure of scaling. The solid line was obtained from the 
recent partial-wave analysis of the single pion photoproduc- 
tion data up to — 2 GeV |6tj)l . while the dashed line was 
obtained from the MAID analysis up to = 1.25 GeV |g3 



C. Possible Substructure of Scaling 

As shown in Figure LT7] 20], though the data with cur- 
rent total uncertainties are consistent with the global 
scaling behavior, the data suggest some substructure for 



both processes: the point around 3.0 GeV is higher than 
those at 2.7 GeV and 3.4 GeV. This might be a hint of 
possible oscillations around the scaling value, similar to 
what was observed in the pp elastic scattering data . 
Future measurements planned at JLab [issj. with finer 
binning in beam energy, are essential for the confirma- 
tion of such oscillatory scaling behavior. 

The substructure or the possible oscillatory scaling be- 
havior in pion photoproduction may arise from the same 
mechanism as in the case of pp scattering, i.e. the inter- 
ference between the long-distance (due to Landshoff di- 
agrams) amplitude and short-distance amplitude |37| , or 
the interference between resonances with a pQCD back- 
ground |38j . But it may also be due to other mechanisms, 
such as high energy resonances around 3 GeV [6l| . the 
interference between the amplitudes associated with dif- 
ferent helicity changes based on the generalized counting 
rule [3ij|4(J,|4l|], or the breaking of the local quark-hadron 
duality above the resonance region [4*2^ . 

The generalized constituent counting rule [3^ was de- 
rived for hard exclusive processes involving parton or- 
bital momentum and hadron helicity flip. It predicts that 
Af( 7 AT T -> ttJVJ - s~ 3 and da/dt^N^ -> ttATJ - s~ 8 
for the helicity- flip case, while M{^N^ — ► -kN-\) ~ s~ 5 / 2 
and da / 'dt(-f N 'j — > nNf) ~ s" 7 for the helicity conserving 
case. By including the amplitudes with helicity changes 
and their interference with the long-distance amplitudes, 
one is able to reproduce |4l| the anomalous oscillations 
of the differential cross section da/dt around the scaling 
value in pp elastic scattering |l5l especially in the low 
energy region (s < 10 GeV 2 ), as well as the very large 
spin-spin correlation [l4| . This may be also true for pion 
photoproduction processes. 

The locality of quark-hadron duality means that the 
local averages of physical variables measured in the res- 
onance region are equal to those measured in the deep- 
inelastic or scaling region. The locality of quark-hadron 
duality can be realized in a simple model of a compos- 
ite system with two spinless charged constituents de- 
scribed by harmonic oscillator wave functions with prin- 
cipal quantum number N and orbital angular momentum 
L (< N). The destructive interference between the high 
density of overlapping resonances leads to the smooth 
scaling behavior at high energies. But for medium en- 
ergies, the locality of quark-hadron duality may break 
down and a sizable oscillation around the scaling value 
may appear above the resonance region due to the orbital 
angular momentum dependence of the resonances [42j . In 
this mechanism, the energy increase reduces the oscilla- 
tion amplitude, and the Q 2 dependence may be nontriv- 
ial. If a subset of resonances are relatively suppressed at 
large Q 2 , there will be significant shifts in the position 
and magnitude of oscillations. The deviation pattern pro- 
duced by the resonance degeneracy breaking requires no 
simple periodicity. The experimental data can thus dis- 
tinguish this mechanism from others. 
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D. Cross Section Enhancement around 2.2 GeV 

Another interesting feature of the differential cross sec- 
tion is an apparent enhancement with a sharp drop for 
both channels of the charged pion photoproduction at 
6 cm. = 90°, at a center-of-mass energy ranging approxi- 
mately from 1.8 to 2.5 GeV, as shown in Figure IT7I Note 
that it is the scaled cross sections s"' da/dt that are plot- 
ted. The drop is even faster in terms of non-scaled cross 
section. A similar cross section enhancement was also 
observed in neutral pion photoproduction |6j| • But 
the enhancement patterns are different for the kinematics 
at 9 c . m . = 70° and 8 c . m . = 50°. 

Some speculation can be made about the enhance- 
ment. It might be due to the known baryon resonances 
around this energy, for example Gi7(2190), Hn(2220) 
and Gig(2250) j6j, just as is the case at lower energies. 
It might relate to some missing resonances |69|. which 
were predicted by the constituent quark model but have 
not been seen experimentally. The value of center-of- 
mass energy hints that the enhancement might be associ- 
ated with the strangeness production threshold, which is 
around 2 GeV to produce a (j> meson of mass 1 GeV. The 
resonances at strangeness and charm production thresh- 
olds were assumed in an approach to explain the strong 
spin-spin correlation and oscillatory scaling in elastic pp 
scattering j3^| . It is worthwhile to mention that a broad 
bump near 2.2 GeV appears in the n~p total cross sec- 
tion, while it is not clear in ir~p elastic cross section |6ll . 
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FIG. 18: (Color online). Exclusive charged pion ratio 



versus center-of-mass energy ^/s and momen- 



d(T / dti^TL — >7T p) 

dtr I dt (7p — >7r+ n) 

turn transfer squre \t\ at pion center-of-mass angle d c .m. ~ 90° 
from JLab E94-104 and previous world data |4g|. together 
with the SAID |M, MAID |67| | and one-hard-gluon-exchange 
calculation [3E5- 



E. Exclusive Charged Pion Ratio 

One can form the exclusive charged pion ratio 
d f / /t ( r i ^+ P \ based on the E94-104 data. As shown 

da / dt(~fp— ^Tr+n) 

in Figure IT8l the exclusive charged pion ratio has some 
oscillations at low energies due to the isospin depen- 
dence of the resonances, which can be described by the 
SAID 66] and MAID calculations available at low en- 
ergies. The big jump around 2 GeV might be associated 
with the isospin-dependent resonances nearby, or with 
the strangeness production threshold (around 2 GeV for 
(j> production). The lowest order (leading- twist) calcula- 
tion based on one-hard-gluon-exchange diagrams [44ll45| . 
which is only valid at high energies, predicts a smooth 

and simple behavior of ffig^g * (^±f^) 2 af- 
ter the nonperturbative components represented by the 
form factors are divided out in the ratio. The theoretical 
prediction seems to agree with the 9 c . m , = 90° data at 
the two highest energies. 

The exclusive charged pion ratio was also calculated to 
higher orders (twist-2 and twist-3) within the handbag 
mechanism considering both quark helicity flip and non- 
flip [jjj. The more precise approach led to the same re- 



|Of| 



suit as the leading- twist prediction when \C%\ » 
The invariant functions \Cf \ (i = 1,4) are the coeffi- 
cients for the four gauge invariant covariants into which 
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FIG. 19: (Color online). Exclusive charged pion ratio 
do-/dt(-fn^Tv p) vergus C enter-of-mass energy Js at different 
pion center-of-mass angles from JLab E94-104. The solid 
curve is calculated by using Equation 2] while the dashed 
one considers the nucleon mass by using Equation 
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the meson photoproduction amplitudes can be decom- 
posed. Both |Cf | and \C%\ contribute only to the quark 
helicity conserving amplitudes while |Cf | and \C[\ gen- 
erate quark helicity flips. The charged pion ratio will 
become infinity at large s if \C^\ is dominant, which is 
clearly not supported by our data. 

The calculation was for massless particle. The effects 
due to nucleon mass may be important at an energy scale 
of few GeV, as estimated by the difference between the 
solid and dashed curves in Figure ED The calculation 
shown by the dashed curve considers the nucleon mass by 
using the identification [jlj of the Mandelstam variables 
with the experimental ones (s oxp , t cxp , u exp ) 

S <S CX p THpi t — ^cxpi — ^cxp (9) 

so that s + t + u ~ 0. FigureH^lalso shows that the agree- 
ment at forward angles is not as good as the case at 90°, 
which might be related to the relatively lower momentum 
transfer. Measurements at higher energies are necessary 
to fully check the theoretical approach, though the cur- 
rent comparison between prediction and data seems to 
suggest that the handbag mechanism is at work for the 
pion photoproduction processes with dominant quark he- 
licity non-flip amplitudes and |Cf | >> |Cf |. 



VI. OUTLOOK 

As mentioned above, a new experiment |r38| was pro- 
posed in JLab Hall A to measure the jn — > ir~p pro- 
cess with deuterium and carbon targets, as well as the 
7P — > TT + n process with a hydrogen target. With very 
fine steps in center-of-mass energy, approximately 0.07 
GeV, the new experiment will be able to elucidate more 
details about the possible substructure of the scaling be- 
havior. In addition, the nuclear transparency of carbon 



in the pion photoproduction process will be measured 
for the first time. This should enable us to test some 
theoretical predictions such as the nuclear filtering effect 
and color transparency. The latter was suggested by the 
helium transparency measurement in E94-104 21]. 

Since the cross section decreases relatively slowly as 
energy increases, the measurements for single pion pho- 
toproduction processes can be greatly extended with the 
JLab 12 GeV upgrade. Since the charm production 
threshold will be crossed, one would be able to inves- 
tigate the resonance at this threshold that was assumed 
in an approach to explain the anomalies in pp scatter- 
ing [3^]. In addition, one would be able to further test 
the scaling and charged pion ratio predictions, especially 
at forward angles. 
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